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Subject of the Research Work. Relativistic nucleus-nucleus collisions are very well suited 
for investigation of excited nuclear matter properties under extreme conditions of high den- 
sities and high temperatures, which are the subject of intense studies both experimentally 
and theoretically. A new collective phenomena have been revealed in high energy heavy ion 
collisions, which allow to study the nuclear equation of state (EOS) and offers an unique 
information about the possible phase transition in the nuclear matter as well. 

The experimental discovery of such states is impossible without understanding the mech- 
anism of collisions and studying the characteristics of multiparticle production in nucleus- 
nucleus interactions. During last 10-15 years several theoretical models have been proposed 
for the study of nucleus-nucleus collisions at high energy. It is assumed that, the intranu- 
clear cascade model (CAScade Intranuclear Model Made in Rossendorf - CASIMIR) and 
Quark Gluon String Model (QGSM) allow to investigate hadron-nuclear, hadron-nucleus 
and nucleus-nucleus collisions at wide energy interval and especially in the range of inter- 
mediate energy y/s < 4 GeV). 

In order to study the dynamics of the relativistic nucleus-nucleus collisions in the se- 
ries of experiments the inclusive spectra of kinematical variables have been obtained for 
secondary protons and mesons. It is interesting to study the correlation between Lorenz 
invariant variables, such as the rapidity Y and transverse momentum P T of the produced 
particles. Interesting information about the collision process and its space-time evolution 
can be obtained from the analysis of secondary particles angular distributions with the aim 
to determine the degree of anisotropy in particles emission. It is also important to investigate 
the energetic spectra of secondary particles at different intervals of emission angle 6>/ ab . 

In the energy range of (0.5 - 10) GeV/nucleon the nuclear matter is compressed and 
heated more than at lower beam energies. Two different signatures of collective flow have 
been predicted by nuclear shock wave models and ideal fluid dynamics: 1) the bounce-off 
of the compressed matter in the reaction plane and 2) the squeeze-out of the participant 
matter out of the reaction plane. Collective flow has been already observed for protons, light 
nuclei, pions, kaons and (A-hyperons in nucleus-nucleus collisions at the energies: of (0.1-1.8) 
GeV/nucleon at BEVALAC (Lawrence Berkeley Laboratory - USA), GSI/SIS (Darmstadt 

- Germany), 3.7 GeV/nucleon at Dubna (Joint Institute Nuclear Research - Russian), 2 

- 14 GeV/nucleon AGS BNL (Brookhaven National Laboratory - USA) and 60 and 200 
GeV/nicleon at CERN SPS (European Organization for Nuclear Research - Switzerland), 
and more recently by the STAR and PHENIX at RHIC BNL at v 7 ^^ =130 CeV/nucleon. 
Many different methods have been proposed for the study of the flow effects in relativistic 
nuclear collisions, of which the most commonly used are the directed transverse momentum 
analysis technique developed by P. Danielewicz and G. Odyniec and the method of the 
Fourier expansion of azimuthal distributions proposed by Demoulins et al. and Voloshin 
and Zang. 

The investigation of the collective effects can be also studied in multiparticle azimuthal 
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correlations, because it gives the information about space-time evolution of the colliding 
system. In order to study these correlations between the protons and between the 7r~ mesons, 
the phase space can be divided into two hemispheres. The division can be performed by use 
of the rapidity and by the emission angle - the forward and backward emitted particles 
and by energy - the slow and fast particles as well. The azimuthal correlations have been 
already observed between protons and between pions in nucleon-nucleus and nucleus-nucleus 
collisions at energies of (0.4-1.8) GeV/nucleon at BEVALAC, 3.7 GeV/nucleon at Dubna 
and 60 and 200 GeV/nucleon at CERN SPS. 
Aims of the Thesis . 

• The study of the properties of 7r~ mesons produced in Mg-Mg collisions at energy of E=3.7 
GeV/nucl at a SKM-200-GIBS set-up of JINR/DUBNA. Investigations of the correlation 
between Lorentz invariant variables: the rapidity Y and the transverse momentum P?. The 
analysis of angular distribution of n~ mesons and the study of the dependence of n~ mesons 
emission anisotropy coefficient a on the kinetic energy. Investigation of emitted n~ mesons 
energetic spectra in different intervals of emission angle 9i ab in the laboratory and centre of 
mass systems (c.m.s.). 

• Investigation of collective flow effects of ir~ mesons produced in central C-Ne, Mg-Mg, C-Cu 
and O-Pb collisions at energy of E=3.7 GeV/nucl by use of transverse momentum technique 
of R Danielewicz and G. Odyniec. The study of the dependence of the flow parameter F, 
as the slope at midrapidity, on the target mass numbers (At)- 

• Investigation of the collective effects by the multiparticle azimuthal correlations between 
protons and between pions in central C-Ne, Mg-Mg, C-Cu and O-Pb collisions. The study of 
the dependence of the correlation function on the angles Aip between the sumed transverse 
vectors for particles emitted in the forward and backward hemispheres. 

• Comparision of the obtained experimental results with predictions of intranuclear cascade 
model (CASIMIR) and Quark Gluon String Model (QGSM). 

Scientific Novelty. 

• A detailed study of pion production in central Mg-Mg collisions have been carried out: 
a) Correlations between Pt and Y have been studied; b) The anisotropy coefficient a of 
7r~ mesons has been obtained from the angular distributions of pions. Its dependence on 
Ekin/E™°f in the c.m.s. has been studied, c) The dependence of the energy spectra of 7r~ 
mesons emitted in an angular interval of 0° < 0i a b < 180° in the laboratory system and 
@c.m.s. — 90° ± 10°, 90° ± 20° emission angles in c.m.s. on Ap, At have been investigated for 
Mg-Mg, He(Li, C), C-Ne, C-Cu, (C,0)Pb collisions. The values of fireball model parame- 
ters, the temperature E^ and velocity (3 of the cluster formed during the collisions, which is 
emitting tt~ mesons have been defined. The obtained results have been compared with the 
predictions of CASIMIR. 

• Using transverse momentum analysis method of P. Danielewicz and G. Odyniec for the 
first time the reaction plane has been defined by pions only and the collective flow effects of 
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7r~ mesons produced in central Mg-Mg collisions have been performed. The similar investi- 
gation for pions produced in central C-Ne, C-Cu and O-Pb collisions has been carried out. 
The dependence of flow parameter F on A^-target mass numbers has been obtained. The 
obtained results have been compared with the predictions of QGSM. 

• The investigation of the multiparticle azimuthal correlations between protons and between 
7r~ mesons in central C-Ne, Mg-Mg, C-Cu and O-Pb collisions has been carried out. The 
dependence of C(Atp) - correlation function between particle groups emitted in the forward 
and backward hemispheres has been studied for the first time. By fitting the experimental 
spectra the asymmetry coefficient £ and the strength of correlation ( has been determined. 
The type of correlations between protons (C-Ne, C-Cu) and between 7r~ mesons (C-Ne, Mg- 
Mg, C-Cu, O-Pb) and the dependence of the correlation parameters £ and ( on the A^-target 
mass numbers have been investigated. The obtained results have been compared with the 
predictions of the QGSM. 

Practical Value . 

• The anisotropy coefficient a of n~ mesons produced in central Mg-Mg collisions at energy of 
3.7 GeV/nucleon increases linearly with E kin / E™£ x starting from approximately 150 MeV. 
In the range up to 150 MeV the coefficient is small a = 0.085 ± 0.020, therefore the ir~ 
mesons are emitted almost isotropically. 

• The shape of the energetic spectra of ix~ mesons emitted at fixed angles in He(Li,C), C- 
Ne, Mg-Mg, C-Cu, and (C,0)Pb collisions in both laboratory (10° < 6 lab . < 180°) and CM. 
(Oc.m.s. = 90° ± 10°, 90° ± 20°) systems does not depend on and A T for all colliding pairs 
of nuclei. The values of the slope parameter E Q of the spectra decrease with increasing of 
9iab from 10° to 110° and remain practically constant at ~50-60 MeV for angles larger than 
110°. 

• For 7T~ mesons produced in He(Li,C), C-Ne, Mg-Mg, C-Cu, and (C,0)Pb collisions the 
angular dependence of E Q is well described by the fireball model. The temperature and 
the velocity (3 of the cluster emitting ir~ mesons have been obtained from the model. The 
temperature weakly depends on Ap, At, while the velocity does not depend on Ap, At- 

• The obtained collective flow effects and multiparticle azimuthal correlations in central C- 
Ne, Mg-Mg, C-Cu and O-Pb collisions at energy of 3.7 GeV/nucleon, show the persistence 
of these collective phenomena in the range of intermediate energy, i.e. from the Bevalac and 
GSI/SIS up to Dubna, AGS, RHIC and SPS energies. 

• The statistical programs of experimental data processing and analyses have been elaborated 
for the personal computers in High Energy Physics Institute of Tbilisi State University, which 
can be used for the processing and analysis of the data of the other experiments. 

Structure and Basic Contents of the Thesis . The thesis contains an introduction, five 
chapter, conclusions and referred bibliography. It consists of 112 pages, accounts 30 figures 
and 7 tables. The references include 147 items. The experimental data obtained at the 
Laboratory of High Energies JINR by the SKM-200-GIBS Collaboration have been used. 
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Introduction presents the subjects of the thesis, its aims and obtained results. Structure 
of the thesis and its short overview are given. 

Chapter I is devoted to the description of the spectrometer GIBS, which has been con- 
structed at the Laboratory of High Energies, at JINR/DUBNA, and is a modified version of 
the SKM-200 set-up. The Facility consists of a 2m streamer chamber with the volume of a 
2x1x0.6 m 3 , placed in a magnetic field of ~0.8 T (~0.9 T for GIBS) and a triggering system. 
In the first Chapter the main elements of both spectrometers are described and the physical 
parameters of spectrometers have been compared. The streamer chamber was exposed to 
beams of 4 He, 12 C, 16 0, 20 Ne (SKM-200) and 12 C, 19 F, 24 Mg (GIBS) nuclei accelerated in 
the synchrophasotron up to at energy per incident nucleon of 3.7 GeV. 

Solid targets Li, C, Cu and Pb in the form of thin discs with thickness of (0.2 ~ 0.5) 
g/cm 2 (the thickness of Li was 1.5 g/cm 2 ) were mounted within the fiducial volume of the 
SKM-200 chamber. Neon-gas filling of the chamber also served as a nuclear target. C, F and 
Mg solid targets with thickness 0.99 g/cm 2 , 1.34 g/cm 2 and 1.56 g/cm 2 , correspondingly, 
were used for GIBS set-up. Experimental set-up and the logic of the triggering system are 
presented in Fig. 1. The "inelastic" trigger, consisting of two sets of scintillation counters 
mounted upstream (SI - S4) and downstream (S5, S6) the chamber, has selected all inelastic 
interactions of incident nuclei on a target. The "central" trigger was selecting events defined 
as those without charged projectile spectator fragments and spectator neutrons (P/Z > 3 
GeV/c), emitted within a cone of half angle 6 ch = 2.4° or 2.9° and 9 n = 1.8° or 2.8° for 
SKM-200, and 9 ch = 6 n = 2.4° for GIBS set-up. 

The streamer chamber pictures were scanned twice at HEPI TSU and JINR/DUBNA. 
The final results of scanning and measurement of events have been recorded on the Data 
Summar Tapes (DST) and then have been analysed by using the standard statistical pro- 
grams HBOOK and PAW. 

Chapter II is devoted to the description of the intranuclear cascade model (CAScade 
Intranuclear Made in Rossendorf) CASIMIR and the Quark Gluon String Model (QGSM) 
wich have been used for comparison with the experimental results presented in this thesis. In 
CASIMIR model the colliding nuclei are represented as an assembly of nucleon, the positions 
and momenta of which are randomly distributed according to realistic nuclear density and 
the law of a cold free Fermi gas, respectively. The nucleon momenta are distributed according 
to a sharp Fermi sphere with a radius of Pp=270 mV/c. The cascade is subsequently traced 
in time by the Monte Carlo technique. The nucleons are assumed to move freely along 
straight lines between collisions. For any pair of particles the minimum distance of approach 
is 



where <Jtot(y/s) is the total cross section for the pair as a function of CM energy. The 
inelastic scattering of nucleons, which is almost entirely dominated by single- and double- 




(1) 
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pion production in the energy range of our interest, is assumed to proceed via the formation 
of two different resonances A(1232) and iV*(1640). 

In the string models the particles are produced via the quark-gluon strings formations 
and decays. The quark-gluon strings are excited objects consist of quark- ant iquark or quark- 
diquark pairs connected via color field. The string masses are formed as a result of composite 
quarks longitudinal motion. The Monte Carlo generator COLLI is based on the Quark Gluon 
String Model. It is possible to study hadron-nuclear, hadron-nucleus and nucleus-nucleus 
collisions in a wide energy interval and especially in the range of intermediate energy (y/s < 
4 GeV) by this program. The procedure of event generation consists of three steps: 1) the 
definition of configurations of colliding nucleons, 2) production of quark-gluon strings and 3) 
fragmentation of strings (breakup) into observed hadrons. The coordinates of nucleons are 
generated according to a realistic nuclear density. It is assumed that the distance between 
them is greater than 0.8 fm. The maximum nucleon fermi momentum is 

P F = (37r 2 ) 1/3 -/>-p 1/3 M, (2) 

where h=0.197 fm-GeV/c, p is nuclear density. For main NN and 7rN interactions the fol- 
lowing topological quark diagrams were used: binary, "undeveloped" cylindrical, diffractive 
and planar. Binary process makes a main contribution. This reaction predominantly results 
in the production of resonances (for instance, NN — > iVA ++ ), which are the main sources 
of pions. The QGSM simplifies the nuclear effects, in particular, coupling of nucleons inside 
the nucleus is neglected, and the decay of excited recoil nuclear fragments and coalescence 
of nucleons is not included. 

Chapter III is devoted to the detailed study of pion production in central Mg-Mg collisions 
at energy of E=3.7 GeV/nucl. The choice of n~ mesons is due to the fact, that they are 
dominantly produced particles, carry information about the dynamics of collisions and can 
be unambiguously identified. Also, the production of n~ mesons is the predominant process 
at the energies of the Dubna synchrophasotron. In this chapter the possible sources of 
experimental biases and appropriate correction procedures for 7r~ mesons are described. It 
has been shown that the dependence of the average kinematical characteristics, < Pt > and 
< Y >, of 7r~ mesons on the multiplicity n_/Ap differs from that for NN collisions at the 
same energy, which is due to nuclear effects. The correlation between < > and Y (Fig. 
2) has been studied. One can see that < Pt > is smaller in the fragmentation regions of the 
projectile and target. The maximum of < Pt > corresponds to yi\T^=l.lA. The shape of the 
Pt distributions in various intervals of Y are described. In the central region of Y the Pt 
spectra are less steep and extended up to 1.4 GeV/c, in the fragmentation regions the spectra 
are more sloping and extended up to (0.5 -j- 0.8) GeV/c. The shape of Pt dependence is 
similar in the fragmentation regions of projectile and target. 

The degree of anisotropy in pion emission can be more directly determined by studying 
the angular distributions in the c.m.s. of colliding nuclei. The approximation of this distri- 
butions by quadratic form gives for anisotropy coefficient: a exper = 0.62 ± 0.02 and a moc i e i = 
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0.63 ± 0.02 (using the CASIMIR code). It has been shown, that the anisotropy coefficient 
increases linearly with E kin / E™°^ (E kin is the kinetic energy, E^^ - the maximum available 
energy in the c.m.s. of 7r~ mesons) starting from approximately 150 MeV (Fig. 3). In the 
range up to 150 MeV the coefficient is small a = 0.085 ± 0.020, therefore the 7r~ mesons 
are emitted almost isotropically, which may be caused by an increase of the number of NN 
collisions. The anisotropy coefficient increases linearly with E k i n jE^}°^ also for the generated 
data, but the slope is more steep than for the experimental data. 

The energy spectra of ir~ mesons at fixed angles in the laboratory and cm. systems in 
He-Li, He-C, C-Ne, Mg-Mg, C-Cu, C-Pb and O-Pb collisions at energy of E =3.7 GeV/nucl 
have been investigated. The spectra of 7r~ mesons depend exponentially on their kinetic 
energy. The values of the slope parameter E (9i a i,) of the spectra - the result of exponential 
approximation of energy spectra are presented in Table 1. One can see that in the forward 
hemisphere (#/ a b. < 90°) two groups of nucleus pairs can be distinguished: approximately 
symmetric and lighter He(Li, C), C-Ne, Mg-Mg and asymmetric and heavier C-Cu, (C, 
0)Pb. In each group E Q parameters coincide within errors and, between these two groups, 
differ no more than 15%. Concerning tx" mesons from the backward hemisphere (9i a b. > 90°), 
E a parameters coincide within errors for all interactions and there is no dependence on Ap, 
A T . The values of the parameter E Q in the cm. system have the same tendency as in the 
laboratory system. The shape of the energy spectra at angular intervals 9 c . m . s . = 90°±10° and 
90° ± 20° does not depend on A P , A T for all pairs of nuclei. The values of E Q decrease from 
(250-300 MeV to ~60 MeV with increasing 6 lab . > 10° from 10° to 110° (Fig. 4). From the 
dependence of E Q parameter on 9i ab . > 10° for each colliding pairs of nuclei, the temperature 
Eq of the cluster formed in the collision and emitting the particles, and its velosity (3 have 
been obtained by use of fireball model approximation E Q = E^ /(l — (3cos9i ab ). The values 
of E^ are comparatively larger ~75 MeV for approximately symmetric and lighter groups 
of nuclear pairs He(Li, C), C-Ne, Mg-Mg, than for asymmetric and heavier pairs of nuclei 
C-Cu, (C, 0)Pb -69 MeV. 

Chapter IV is devoted to the experimental results obtained from the in-plane transverse 
momentum analysis of 7r~ mesons in central C-Ne, Mg-Mg, C-Cu and O-Pb interactions at 
energy of E=3.7 GeV/nucl. As the ir~ mesons are emitted mainly from decays of A-isobars 
(at least ~ 80%), we decided to investigate whether a single 7r~ meson of the event knows 
something about its origin and the question whether they are collectively correlated. Many 
different methods were proposed for the study of flow effects in relativistic nuclear collisions, 
of which the most commonly used are - the transverse momentum analysis technique devel- 
oped by P. Danielewicz and G. Odyniec and Fourier analysis of the azimuthal distributions 
on the event-by-event basis in relatively narrow rapidity windows proposed by M. Demoulins, 
S. Voloshin and Y. Zhang. We have used the technique of P.Danielewicz and G. Odyniec 
in our analysis. The method involves two basic ideas: 1) to select the rapidity range and 
rapidity dependent waiting factors to the real reaction plane, and 2) to remove trivial and 
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spurious self-correlations from the projections. Autocorrelations are removed by calculating 
Q individually for each particle without including that particle into sum: 

(Tj=T,Ui Pti (3) 

where i is a particle index and Ui is a weight, taken as 1 for y^> j cms and -1 for ji< j cms 
{jams is c.m.s. rapidity and ji is the rapidity of particle i). The transverse momentum of 
each particle in the estimated reaction plane is calculated as: 

P xj '=£wi-(P?rP£)/\$\ ( 4 ) 

The Q vector have been constructed from the Pn of only 7r~ mesons event-by- event for 
the events with multiplicity n_ > 7 in Mg-Mg collisions. It is known, that the estimated 
reaction plane differs from the true one, due to the finite number of particles in each event. 
The component P x in the true reaction plane is systematically larger then P x ' the component 
in the estimated plane, hence < P x >=< P x ' > / < costp >, where is the angle between 
the estimated and true planes and the correction factor K — 1/ < cosip >. According to 
the method of Danielewicz and Odyniec, for the definition of < cosip >, we divided each 
event randomly into two equal sub-events, constructed vectors Qj and Qu and estimated 
azimuthal angle 4>i,2 between these two vectors < cosip >=< cosipi^ >• The correction 
factor, averaged over all the multiplicities, K = 1.51 ± 0.05 has been obtained. Figure 5 
shows the dependence of the estimated < P X '(Y) > on Y for pions produced in the central 
Mg-Mg collisions. The data exhibit S-shape behavior similar to the form of the < P x > 
spectra for protons and pions obtained at lower energies and identified as nucleon and pion 
collective flow. The slope at midrapidity has been extracted from a linear fit of the data 
for Y between 0.2 2, the value of F = 48 ± 5 MeV/c have been obtained. In the model 
calculations the reaction plane is known a priori and is refered as the reaction plane. For 
simulated events (QGSM) the component in the true reaction plane P x has been calculated 
and the dependence of < P X '{Y) >on7 has been obtained for both not fixed b (6225 events) 
and fixed < b >=1.34 fm (6212 events) impact parameters (Fig. 5). The experimental and 
QGSM results coincide within errors. The values of F, obtained from the QGSM are: F = 
53 ± 3 (MeV/c) - for not fixed b, and F = 51 ± 4 MeV/c - for 1.34 fm. 

The similar in-plane analysis has been performed on 7r~ mesons too in central C-Ne, C-Cu 
and O-Pb collisions at the same energy. In this case the reaction plane has been obtained by 
only 7r~ mesons. The observed dependences show the typical S-shape behavior reflecting the 
presence of flow effects. The values of flow F increases with the mass number of target At 
from 33 ± 4 (MeV/c) to 50 ± 6 (MeV/c) (see Table 2). The obtained experimental results 
in central C-Ne, C-Cu and O-Pb collisions have been compared with the predictions of the 
QGSM. For this purpose C-Ne (6272), C-Cu (9327) and O-Pb (2431) collisions have been 
generated for impact parameters b=2.20, 2.75 and 3.74 (fm), respectively. The Quark Gluon 
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String Model reproduces the < P x > distributions and satisfactorily describes the obtained 
experimental results for all pairs of nuclei. 

It seems, that the obtained results indicate that the flow behavior of n~ mesons is due 
both to the flow of A resonances and of the nuclear shadowing effect. Our experimental 
results, obtained using the streamer chamber technique, provide quantitative information 
on the directed flow of 7r _ mesons and their dependence on target-mass, complementing the 
experimental data available from the BEVALAC, GSI-SIS, AGS and CERN/SPS. 

Chapter V is devoted to the multiparticle azimuthal correlations between protons and be- 
tween pions in central C-Ne, Mg-Mg, C-Cu and O-Pb collisions at energy of 3.7 GeV/nucleon. 
Identification of secondary particles becomes very important for the study of products of 
nucleus- nucleus (Ap - At) collisions via collective variables, which are defined event by 
event. The admixture of 7r~ mesons amongst the charged positive particles is about (25 -j- 
27). The statistical method has been used for the identification of 7r + mesons. The main 
assumption is based on the similarity of spectra of n~ and ir + mesons Pt, Pl)- The two- 
dimensional distribution of Pt, Pl variables have been used for identification of n~ mesons. 
After performed identification the admixture of 7r~ mesons amongst the protons is not ex- 
ceeding (5 4- 7)%. The azimuthal correlation function C(Acp) was defined by the relative 
angle between the transverse momentum vectors sums of particles emitted in forward and 
backward hemispheres: 

Qf= J2 Pri (5) 

Vi><y> 

and 

Qb= E Pri (6) 

Vi«y> 

where < y > is the average rapidity in each event. The correlation function C(Atp) is 
constructed as follows: 

C(Aip) = dN/dAv? (7) 

where Aip is the angle between these vectors: 

Aip — arccos(Q^ ■ Qf)/\Qb \ • \Qf\- (8) 

For protons a "back-to-back" negative correlations were observed in C-Ne and C-Cu collisions 
i.e., protons are preferentially emitted at Aip = 180 °. To quantify these experimental results, 
the data were fitted by C(Atp) = 1 +l;cos(Aip), where £ is the asymmetry coefficient. The 
strength of correlation ( was defined as 

C = C(0°)/C(180°) = (i + o/(i-0- (9) 
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In Table 3 the asymmetry coefficient £ < and the strength of correlation £ < 1 for protons 
in both C-Ne and C-Cu interactions are presented. One can see that |£| increases and ( 
decreases with increase of target mass. A back-to-back pion correlations have been obtained 
for a light, symmetric pairs of nuclei C-Ne and Mg-Mg (n_ > 7), where £ and ( parameters 
have the same behavior as for protons (Fig. 6a). For heavy pairs of nuclei C-Ne and O-Pb, 
side-by-side (positive) pion correlations were observed i.e. in this case they are preferentially 
emitted at Atp = 0° (Fig. 6b). The asymmetry coefficient £ and the strength of correlations 
( increase with increase of Ap and At- The variation of pion correlations (from negative for 
light systems of C-Ne, Mg-Mg to positive for heavy systems of nuclei C-Cu, O-Pb), which 
have been observed in our experiment, indicate, that at Dubna energies the flow behavior 
of 7r~ mesons in a light system is due the flow of A resonances (ttN — > A and AiV — > NN), 
whereas the antiflow behavior in a heavier system is the result of the nuclear shadowing 
effect. 

To be convinced that the azimuthal correlations between protons and between pions are 
due to correlations between these particles and can not be the result of detector biases or 
finite- multiplicity effects, we obtained the of C(Aip) on ip for secondaries, where ip is the 
angle between the transverse momentum of each particle emitted in the backward (forward) 
hemisphere and Qf(Qb) vector, respectively. No correlations have been obtained for C-Cu 
interactions both for protons and for pions. Similar results have been obtained for C-Ne 
collisions. 

In order to extend these investigations, the relation between < P x 2 > and the angle Atp, 
where Aip is the opening angle between Q F and Qp vectors has been obtained. One can see 
from Fig. 7, that for protons in C-Ne and C-Cu collisions the distributions show S-shape 
behavior and slopes of distributions increase with target mass; We stress, that at Aip = 90° 
the values of < P x > do not depend on A T . 

The Quark Gluon String Model (QGSM) satisfactorily describes the obtained results for 
protons and for 7r_ mesons for all pairs of nuclei. 

The comparison of our experimental results with the data of various collaborations (BE- 
VALAC, JINR, CERN/SPS) has been carried out and had been showed that these multi- 
particle azimuthal correlations between the same types of secondaries exist for the energy 
interval of (0.1 + 200) GeV/nucleon. 

Conclusions list main results of the conducted research work, which are outlined bellows. 

The Main Results of the Thesis 
1. A detailed study of pion production in central Mg-Mg collisions at energy of E=3.7 GeV 
per nucleon (P = 4.3 GeV/c/nucl) have been carried by the Lorentz invariant variables: the 
rapidity Y and the transverse momentum P?. The average values of < Y > and < P? > 
decrease with increasing of multiplicity, ri-/Ap, different to NN collisions, where the values 
of both average kinematical characteristics remain constant with increase of ri-jAp. n^/Ap 
is the measure of the impact parameter b. 
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The correlation between < P T > and Y has been studied: < Pt > is smaller in the 
fragmentation regions of the projectile and target (the maximum of < Pt > corresponds to 
Dnn = 1.14). In the central region of Y the Pt distributions are less steep and extended up 
to 1.4 GeV/c. 

2. In the analysis of angular distributions of 7r~ mesons in Mg-Mg collision the anisotropy 
coefficient a exper , = 0.62 ± 0.02 has been obtained; The obtained results qualitatively agree 
with the predictions of the intranuclear cascade model CASIMIR a mo d e i = 0.63 ± 0.02. The 
anisotropy coefficient increases linearly with Ek in /E^ x {E^in is the kinetic energy, E^ x - 
the maximum available energy in the c.m.s. of 7r~~ mesons) starting from approximately 150 
MeV. In the range up to 150 MeV the coefficient is small a = 0.085 ± 0.020, therefore the 7r~ 
mesons are emitted almost isotropically. The anisotropy coefficient increases linearly with 
Ekin/ EklT a l so f° r the generated data (CASIMIR), but the slope is more steep than for the 
experimental data. 

3. The energy spectra of 7r~ mesons emitted at fixed angles in the laboratory (10° < 9\ ah < 
90°) and cm. systems (9 c . m . s . = 90°±10° and 90°±20°) in He-Li, He-C, C-Ne, Mg-Mg, C-Cu, 
C-Pb and O-Pb collisions collisions at energy of E=3.7 GeV/nucl depend exponentially on 
their kinetic energy. The shape of spectra in both systems does not depend on A P) A T for all 
colliding pairs of nuclei. The values of the slope parameter E Q of the spectra depend on A P , 
A T for angles 9i ab < 90° and are independent for 9i ah > 90°. The values of E Q decrease from 
~250-300 MeV to ~60 MeV with increasing 9[ a b from 10° to 110°, and remain practically 
constant at ~50-60 MeV for angles larger 110°. The angular dependence of E Q parameter 
has been approximated by the fireball model. From the approximation the parameters of 
the model - the temperature and velosity f3 of the cluster emitting n~ mesons have been 
extracted. The temperature weakly depends on A P , A T while the velocity does not depend 
on Ap, At- 

4. The investigation of in-plane flow effects for n~ mesons in central Mg-Mg collisions by 
the transverse momentum technique of P. Danielewicz and G. Odyniec has been carried out. 
For the first time, the reaction plane heve been determined only by 7r~ mesons in each event. 
The similar analysis of in-plane flow effects for 7r~ mesons production in central C-Ne, C-Cu 
and O-Pb collisions at the same energy has been performed. The observed dependence of 
< P X '(Y) > on Y shows S-shape behavior. The values of F flow parameter, defined as the 
slope at midrapidity, have been extracted for all colliding pairs of nuclei. F increases with 
the mass numbers of target At from 33 ± 4 (MeV/c) to 50 ± 6 (MeV/c). The Quark Gluon 
String Model (QGSM) reproduces the < P x > distributions and satisfactorily describes 
obtained experimental results for all pairs of nuclei. 

5. The multiparticle azimuthal correlations between protons and between pions in central C- 
Ne, Mg-Mg, C-Cu and O-Pb collisions have been studied. The dependence of the azimuthal 
correlation function C(Atp) on the A(p, where Aip, is angle between the transverse momentum 
vectors sums of particles emitted in forward and backward hemispheres, in these collisions 
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have been investigated: 

For protons a "back-to- back" (negative) correlations were observed in central C-Ne and 
C-Cu interactions. The asymmetry coefficient |£ (£ < 0) increases and the strength of 
correlation ( (( < 1) decreases with increase of the target mass At- 

"Back-to-back" pion correlations have been obtained for a light pairs of nuclei (C-Ne and 
Mg-Mg), where £ and ( parameters have the same behavior as for protons. For heavy pairs 
of nuclei (O-Pb and C-Cu) " side- by-side" (positive) pions correlations were observed. The 
asymmetry (£ > 0) and strength of correlation (( > 1) increase with increase of Ap and At- 

The QGSM satisfactorily describes multiparticle azimuthal correlations of protons and 
7r~ mesons for all pairs of nuclei. 

Approbation of the Work. The basic results of the investigations were reported at: Semi- 
nars in High Energy Physics Institute of Tbilisi State University and in Dubna Joint Institute 
of Nuclear Research, International Symposium in Elemantary Particle Physics dediceated to 
the memory G. Chikovani, 21-23 September, 1998, Tbilisi; On the 1999 European School of 
young scientist in High-Energy Physics, Casta-Papiernicka, Slovac Republic, 22 August - 4 
September and on the Seventh International Conference in Nucleus-Nucleus Collisions 2000, 
July 3-7, Strasbourg, France. 
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Table 1. The values of the paramether E - the result of approximating 7r_ mesons energy 
spectra by formula o inv _ = Aexp(E fein /E ). 



ft, , 

Vlab 

i n't" c*v\rz\ 1 q 
lllliCl Velio 

(grad) 


^ ^kin 

I VIA V 1 


E (GeV) 


He(Li,C) 


C-Ne 


Mg-Mg 


C-Cu 


(C,0)Pb 


10-20 


0.2-0.3 


0.408 ± 0.009 


0.390 ± 0.010 


0.378 ± 0.006 


0.349 ± 0.008 


0.318 ± 0.008 


20-30 


0.15-2.1 


0.289 ± 0.007 


0.283 ± 0.009 


0.271 ± 0.005 


0.259 ± 0.006 


0.234 ± 0.007 


30-40 


0.10-1.6 


0.211 ± 0.006 


0.200 ± 0.007 


0.210 ± 0.005 


0.182 ± 0.004 


0.175 ± 0.005 


40-50 


0.05-1.2 


0.161 ± 0.005 


0.159 ± 0.006 


0.170 ± 0.004 


0.140 ± 0.004 


0.136 ± 0.004 


50-70 


0.05-1.0 


0.123 ± 0.004 


0.121 ± 0.006 


0.122 ± 0.003 


0.109 ± 0.003 


0.098 ± 0.003 


70-90 


0.05-0.7 


0.083 ± 0.005 


0.080 ± 0.005 


0.089 ± 0.004 


0.079 ± 0.003 


0.071 ± 0.003 


90-110 


0.05-0.7 


0.059 ± 0.005 


0.057 ± 0.005 


0.059 ± 0.002 


0.059 ± 0.002 


0.058 ± 0.002 


110-130 


0.05-0.5 


0.064 ± 0.011 


0.060 ± 0.008 


0.053 ± 0.005 


0.052 ± 0.003 


0.055 ± 0.003 


130-180 


0.05-0.5 


0.046 ± 0.004 


0.052 ± 0.007 


0.052 ± 0.004 


0.046 ± 0.003 


0.041 ± 0.003 
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Table 2. The number of experimental events N event , the average multiplicity of protons 
(participant) < n prot > and tt~ mesons < n n _ >, and the flow parameter in central C-Ne, 
Mg-Mg, C-Cu and O-Pb collisions (the reaction plane has been obtained by 7r~ mesons 
only). 



A P - A T 


C-Ne 


Mg-Mg 


C-Cu 


O-Pb 


Nevent 


723 


6239 


1866 


732 




12.4 ± 0.5 




19.5 ± 0.6 




< n 7T _ > 


7.8 ± 0.4 


8.2 ± 0.4 


6.6 ± 0.3 


9.8 ± 0.5 


F (MeV/c) 


33 ± 4 


48 ± 5 


41 ± 5 


50 ± 6 
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Table 3. The number of experimental events iV even t and of participant protons iV prot ., 
the asymmetry coefficient (£), the strength of the correlation (() and the average rapidity 
(< y >) of protons in central C-Ne and C-Cu collisions. 



A P - A T 


^ event 






C 


< y > 


C-Ne 


723 


9201 


-0.23 ± 0.05 


0.63 ± 0.09 


1.07 ± 0.07 


C-Cu 


663 


12715 


-0.35 ± 0.05 


0.48 ± 0.06 


0.73 ± 0.05 
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GrIBS 




"INELASTIC" TRIGGER = s i * S 2 * $3* \ * s 5 * % 
"CENTRAL" TRIGGER = y V S 3^ V V^ch 

Figure 1: Experimental set-up. The trigger and trigger distances are not to scale. 
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Figure 2: The dependence of the average transverse momentum < P T > 
on the rapidity Y for n~ mesoms. 
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Figure 3: The dependence of the anisotropy coefficient a on the 
Ekm/E^ x (Ekin is the kinetic energy, E^^ - maximum available en- 
ergy in the c.m.s. of 7r~ mesons): o - experiment, A - CASIMIR. The 
lines are the result of a linear approximation of the experimental and 
generated data, respectively. 
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Figure 4: The angular dependence of the parameter E for A - He(Li, 
C), * - C-Ne, • - Mg-Mg, o - C-Cu, o - (C,0)Pb collisions (P=4.5 
GeV/c/nucl) and * - He-C, + - C-C collisions (P=4.2 GeV/c/nucl). The 
solid curve is the result of approximating the averaged values E{6i ab ) by 
function E Q = / (1- (3 cos 6i ab ) 



20 




Figure 5: The dependence of < P X '{Y) > on Yi ab . for n~ mesons produced 
in the central Mg-Mg collisions: o - the experimental data and QGSM 
generated data for both A - fixed b= 1.34 fm and * - not fixed b impact 
parameters. The solid line is the result of the the linear approximation 
of experimental data in the interval of Y - 0.2 2.0. The curves for 
visual presentation of the QGSM events (solid - for fixed b, dashed -for 
b) - result of approximation by 4-th order polynomial. 
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Figure 6: The correlation function C(A<p) as descreibed in the text for ir -mesons from: 
(a) o - C-Ne, • - Mg-Mg (n_ > 7), o - Mg-Mg (n_ < 7) and (b) A - C-Cu, © - O-Pb 
collisions, x, * - the QGSM generated data, respectively. The solid and dashed curves - 
results of the approximation of the data by function C(A<p)=l + £cos(A<y?). 
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Figure 7: The dependence of < P x > 2 on the Acp for protons from: o - C- 
Ne and A - C-Cu collisions and x, * - the QGSM data, respectively. The 
dashed lines are the result of the linear approximation of experimental 
data in the interval of Aip 60 145 and 45 -h 155, respectively. The 
solid curves for visual presentation - result of approximation by 4-th 
order polynomial. 
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